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Using the formahsm of spin-weighted functions we present an all-sky analysis of polarization 
in the Cosmic Microwave Background (CMB). Linear polarization is a second-rank symmetric and 
traceless tensor, which can be decomposed on a sphere into spin ±2 spherical harmonics. These are 
the analog of the spherical harmonics used in the temperature maps and obey the same completeness 
and orthogonality relations. We show that there exist two linear combinations of spin ±2 multipole 
. . . moments which have opposite parities and can be used to fully characterize the statistical properties 

' of polarization in the CMB. Magnetic- type parity combination does not receieve contributions from 

scalar modes and does not cross-correlate with either temperature or electric-type parity combina- 
tion, so there are four different power spectra that fully characterize statistical properties of CMB. 
We present their explicit expressions for scalar and tensor modes in the form of line of sight integral 
JH , solution and numerically evaluate them for a representative set of models. These general solutions 

differ from the expressions obtained previously in the small scale limit both for scalar and tensor 
modes. A method to generate and analyze all sky maps of temperature and polarization is given and 
the optimal estimators for various power spectra and their corresponding variances are discussed. 
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I. INTRODUCTION 



The field of CMB anisotropies has become one of the main testing grounds for the theories of structure formation and 
early universe. Since the first detection by COBE satellite ^ there have been several new detections on smaller angular 
scales (see for a recent review). There is hope that future experiments such as MAP |^ and COBRAS/SAMBA 0| 
will accurately measure the anisotropies over the whole sky with a fraction of a degree angular resolution, which will 
help to determine several cosmological parameters with an unprecedented accuracy . Not all of the cosmological 
Q ■ parameters can be accurately determined by the CMB temperature measurements. On large angular scales cosmic 
^ \ variance (finite number of multipole moments on the sky) limits our ability to extract useful information from the 
^ ■ observational data. If a certain parameter only shows its signature on large angular scales then the accuracy with 
which it can be determined is limited. For example, contribution from primordial gravity waves, if present, will only be 
important on large angular scales. Because both scalar and tensor modes contribute to the temperature anisotropy one 
cannot accurately separate them if only a small number of independent realizations (multipoles) contain a significant 
contribution from tensor modes. Similarly, reionization tends to uniformly suppress the temperature anisotropies for 
all but the lowest multipole moments and is thus almost degenerate with the amplitude [^,^ . 

It is clear from previous discussion that additional information will be needed to constrain some of the cosmological 
parameters. While the epoch of reionization could in principle be determined through the high redshift observations, 
primordial gravity waves can only be detected at present from CMB observations. It has been long recognized that 
there is additional information present in the CMB data in the form of linear polarization . Polarization could be 

particularly useful for constraining the epoch and degree of reionization because the amplitude is significantly increased 
and has a characteristic signature Recently it was also shown that density perturbations (scalar modes) do not 
contribute to polarization for a certain combination of Stokes parameters, in contrast with the primordial gravity 
waves which can therefore in principle be detected even for very small amplitudes. Polarization information 

which will potentially become available with the next generation of experiments will thus provide significant additional 
information that will help to constrain the underlying cosmological model. 

Previous work on polarization has been restricted to the small scale limit (e.g. |8-10 r^ , p^ , p^ ) . The correlation 



functions and corresponding power spectra were calculated for the Stokes Q and U parameters, which are defined 
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with respect to a fixed coordinate system in the sky. While such a coordinate system is well defined over a small 
patch in the sky, it becomes ambiguous once the whole sky is considered because one cannot define a rotationally 
invariant orthogonal basis on a sphere. Note that this is not problematic if one is only considering cross-correlation 
function between polarization and temperature p] , ^ , where one can fix Q or ?7 at a given point and average 
over temperature, which is rotationally invariant. However, if one wants to analyze the auto-correlation function of 
polarization or perform directly the power spectrum analysis on the data (which, as argued in fl^ , is more efficient 
in terms of extracting the signal from the data) then a more general analysis of polarization is required. A related 
problem is the calculation of rotationally invariant power spectrum. Although it is relatively simple to calculate Q and 
U in the coordinate system where the wavevector describing the perturbation is aligned with the z axis, superposition 
of the different modes becomes complicated because Q and U have to be rotated to a common frame before the 
superposition can be done. Only in the small scale limit can this rotation be simply expressed |Q, so that the 
power spectra can be calculated. However, as argued above, this is not the regime where polarization can make most 
significant impact in breaking the parameter degeneracies caused by cosmic variance. A more general method that 
would allow to analyze polarization over the whole sky has been lacking so-far. 

In this paper we present a complete all-sky analysis of polarization and its corresponding power spectra. In section 
§2 we expand polarization in the sky in spin- weighted harmonics | [r^ , pO| , which form a complete and orthonormal 
system of tensor functions on the sphere. Recently, an alternative expansion in tensor harmonics has been presented 
pGf . Our approach differs both in the way we expand polarization on a sphere and in the way we solve for the 
theoretical power spectra. We use the line of sight integral solution of the photon Boltzmann equation |2^ to obtain 
the correct expressions for the polarization-polarization and temperature-polarization power spectra both for scalar 
(§3) and tensor (§4) modes. In contrast with previous work the expressions presented here are valid for any angular 
scale and in §5 we show how they reduce to the corresponding small scale expressions. In section §6 we discuss how 
to generate and analyze all-sky maps of polarization and what is the accuracy with which one can reconstruct the 
various power spectra when cosmic variance and noise are included. This is followed by discussion and conclusions in 
§7. For completeness we review in Appendix the basic properties of spin-weighted functions. All the calculations in 
this paper are restricted to a flat geometry. 



II. STOKES PARAMETERS AND SPIN-S SPHERICAL HARMONICS 



The CMB radiation field is characterized by a 2 x 2 intensity tensor /y . The Stokes parameters Q and U are defined 
as Q = (111 — /22)/4 and U = /12/2, while the temperature anisotropy is given by T = {In + /22)/4. In principle 
the fourth Stokes parameter V that describes circular polarization would also be needed, but in cosmology it can 
be ignored because it cannot be generated through Thomson scattering. While the temperature is invariant under a 
right handed rotation in the plane perpendicular to direction n, Q and U transform under rotation by an angle ip as 

Q' = Qcos2V' + Usm2ip 

U' ^-Qs\tl2iI) + Ucos2iP (1) 

where ei' = cosV-'e'i -I- smipe-^ and 62' = — sin -061 -I- cos?/'e2. This means we can construct two quantities from the 
Stokes Q and U parameters that have a definite value of spin (see Appendix for a review of spin-weighted functions 
and their properties), 

(Q ± iU)'{h) - eT2,:v.(g ± iU){n). (2) 
We may therefore expand each of the quantities in the appropriate spin-weighted basis 

T{n) = ^aTjmYhn(n) 

Im 

2Yim{n) 

Im 

{Q - iU){h) = ^ a-2.im -2yim{rt). (3) 

Im 

Q and U are defined at a given direction n with respect to the spherical coordinate system {eg,e^). Using the first 
equation in ( |A5D one can show that the expansion coefficients for the polarization variables satisfy alg — o-2,i-m- 
For temperature the relation is im ~ ^T.i-m- 
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The main difficulty when computing the power spectrum of polarization in the past originated in the fact that 
the Stokes parameters are not invariant under rotations in the plane perpendicular to h. While Q and U are easily 
calculated in a coordinate system where the wavevector k is parallel to z, the superposition of the different modes 
is complicated by the behaviour of Q and U under rotations (equation For each wavevector k and direction on 
the sky n one has to rotate the Q and U parameters from the k and n dependent basis into a fixed basis on the sky. 
Only in the small scale limit is this process well defined, which is why this approximation has always been assumed 
in previous work [p|-p^ 14 1^. However, one can use the spin raising and lowering operators d and 3 defined in 
Appendix to obtain spin zero quantities. These have the advantage of being rotationally invariant like the temperature 
and no ambiguities connected with the rotation of coordinate system arise. Acting twice with d , d on Q zL ill in 
equation (|) leads to 



3\Q + iU)in)^J2 

Im 



(/ + 2)! 



{1-2)1 
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1/2 
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The expressions for the expansion coefficients are 



aT,im 

a,2,lm 



dnYC^{n)T{n) 
dn2YiUn){Q + iU){n) 

dn YCM3\Q + iU)in) 



a-2,im — 



{1 + 2)IY^/^ 
a~2)!. 
dn -2Y;„,{n){Q - tU){n) 



{1 + 2)1 



-1/2 



dn Yr^{h)d\Q-iU){h). 



Xl-2)\_ 

Instead of 02, im, o,-2,im it is convenient to introduce their linear combinations |2C 

aE,l7n — ~{a2,lm + a-2,im)/2 

aB,lm = i{a2,im — a-2Jm)/2. 



(5) 



(6) 



These two combinations behave differently under parity transformation: while E remains unchanged B changes the 
sign [ pO| , in analogy with electric and magnetic fields. The sign convention in equation makes these expressions 
consistent with those defined previously in the small scale limit . 

To characterize the statistics of the CMB perturbations only four power spectra are needed, those for T, E, B and 
the cross correlation between T and E. The cross correlation between B and E ot B and T vanishes because B has 
the opposite parity of T and E. We will show this explicitly for scalar and tensor modes in the following sections. 
The power spectra are defined as the rotationally invariant quantities 
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in terms of which. 
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a-T.ln 



{^*B,l'm''^BAm) = CBlSl'lSm'm 
{a*B,l'm'^E,lm) = {a*B J' m' ,lm) = 0' 



(8) 



For real space calculations it is useful to introduce two scalar quantities and B{fi) defined as 



E{h) = [3\Q + iU)+d\Q-iU) 



E 

Irn 



{l + 2)\ 



1/2 



B{n) 



Irn 



(/-2)!J °'E,lmYiyn{fl) 

\6^{Q + iU)-6^{Q-iU) 



(9) 



These variables have the advantage of being rotationally invariant and easy to calculate in real space. These are not 
rotationally invariant versions of Q and U, because d ^ and d ^ are differential operators and are more closely related 
to the rotationally invariant Laplacian of Q and U . In I space the two are simply related as 



{i + 2y. 
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0'{E,B),lm- 
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III. POWER SPECTRUM OF SCALAR MODES 



The usual starting point for solving the radiation transfer is the Boltzmann equation. We will expand the perturba- 
tions in Fourier modes characterized by wavevector k. For a given Fourier mode we can work in the coordinate system 
where A; || i; and (ei, 62 ) = (e g, e^). For each plane wave the scattering can be described as the transport through a 
plane parallel medium |22,E3] . Because of azimuthal symmetry only Q Stokes parameter is generated in this frame 
and its amplitude only depends on the angle between the photon direction and wavevector, 11 — n ■ k. The Stokes 
parameters for this mode are Q = Ap (r, k, n) and U — 0, where the superscript S denotes scalar modes, while the 

temperature anisotropy is denoted with Aj, {T,k,fi). The Boltzmann equation can be written in the synchronous 
gauge as @,|| 



AAp 



ikfiAi!^'^ = --h - -(h + 6ri)P2(fi) + k 



-A, 



■ ikfiA^p^ 



AS) 



n 



A 



is) 

P2 



[i-^^2(A^)]n 



(11) 



Here the derivatives are taken with respect to the conformal time r. The differential optical depth for Thomson 
scattering is denoted as k = aUgXeCrT^ where a(r) is the expansion factor normalized to unity today, rie is the electron 
density, Xe is the ionization fraction and ax is the Thomson cross section. The total optical depth at time r is obtained 
by integrating k, k{t) — k{T)dT. The sources in these equations involve the multipole moments of temperature 
and polarization, which are defined as A(fc,/z) = '}2ii{2l + 1)(— i)'A;(fc)P;(/i), where Pi{ijl) is the Legendre polynomial 
of order /. Temperature anisotropics have additional sources in metric perturbations h and r/ and in baryon velocity 
term w^. 

To obtain the complete solution we need to evolve the anisotropics until the present epoch and integrate over all 
the Fourier modes. 
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i[/(^))(n) 
iC/(^))(n) 



y d3fce(fc)e-2»^--A(f)(T = To,fc,M) 
y d3fce(fc)e2'*^."A(f)(T = To,fc,M), 



(12) 



where is the angle needed to rotate the k and h dependent basis to a fixed frame in the sky. This rotation 
was a source of comphcations in previous attempts to characterize the CMB polarization. We will avoid it in what 
follows by working with the rotationally invariant quantities. We introduced ^(fc), which is a random variable used 
to characterize the initial amplitude of the mode. It has the following statistical property 



(r(fei)e(fe2)> = P0(fc)<5(fci-fc2), 

where P^{k) is the initial power spectrum. 

To obtain the power spectrum we integrate the Boltzmann equation ([l^) along the line of sight |21 



(13) 



A, 



(To,fc,/^)= rdTe*-'^4^)(fc,T) 



dre"''g(r)n(/c,r) 



4^^(fc,T) =g { AT,o + 2a + 



+ e ''(?) + &) + g \a + — 




n 



— L^rp^ 
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PO ' 



(14) 



where x = k(TQ — r) and a = {h + 6r])/2k^. We have introduced the visibility function g^r) — kexp{—K). Its peak 

defines the epoch of recombination, which gives the dominant contribution to the CMB anisotropies. 

Because in the fc || i: coordinate frame U — and Q is only a function of jj. it follows from equation (A3) that 

d ^{Q + ill) ~ d'^[Q — ill), so that 2a/„i — ^2^/™- Scalar modes thus contribute only to the E combination and B 
vanishes identically. Acting with the spin raising operator twice on the integral solution for A^'' (equation |lj) leads 
to the following expressions for the scalar polarization E 



A^^\ro,k,^i) 



r drg{r)n{T,k)dl[{l-f,'re^^^] 
Jo 

dTg{T)n{T,k) {l + dlf{x\'^^). 



(15) 



The power spectra defined in equation (|7|) are rotationally invariant quantities so they can be calculated in the 
frame where fc || 2 for each Fourier mode and then integrated over all the modes, as different modes are statistically 
independent. The present day amplitude for each mode depends both on its evolution and on its initial amplitude. 
For temperature anisotropy T it is given by [pT| 



a 



is) 

Tl 



1 



21 + 1 



d^kp^{k)j2 



dilYCM I drS^^'\k,T) e"^ 



= (47r)2 / k^dkP^{k) 



dTS!^\k,T)ji (x) 



(16) 



where ji{x) is the spherical Bessel function of order / and we used that in the || 2: frame / dilYi^{n) e^^^ = 
■y/47r(2^ + l)i''ji{x)Smo- For the spectrum of E polarization the calculation is similar. Equation ( p^ ) is used to 
compute the power spectrum of E which combined with equation ([l0|) gives 
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2/ + 1 (^ + 2)! 

,{l-2)\ 



I fj dnY.Ufi) £ drgirMk, r) ([1 + dl^x^e^^n 



(47r)^|±|| / k'dkP,{k) 



3 r° 

3 '■^o 



drg(T)n(r,fc) {[1 + d'^nx'jiix)] 

2 



dr9(r)n(T,A:):^ 



(17) 



To obtain the last expression we used the differential equation satisfied by the spherical Bessel functions, j'/ + 'ij'Jx + 
[1 - /(/ + l)/x'^]ii = 0. If we introduce 



5f (fcr) 



dTS!^'\k,T)ji{x) 



3g(T)n(T,fc) 
4a;2 



then the power spectra for T and -B and their cross-correlation are simply given by 

Ch'^ = {^^f j k'dkP4k)A^^',\k)A'^i\k). 
Equations (^8|) and (|l^) are the main results of this section. 



(18) 



(19) 



IV. POWER SPECTRUM OF TENSOR MODES 



The method of analysis used in previous section for scalar polarization can be used for tensor modes as well. The 
situation is somewhat more complicated here because for each Fourier mode gravity waves have two independent 
polarizations usually denoted with + and x . For our purposes it is convenient to rotate this combination and work 
with the following two linear combinations, 



(20) 



where f 's are independent random variables used to characterize the statistics of the gravity waves. These variables 
have the following statistical properties 



where Ph{k) is the primordial power spectrum of the gravity waves. 

In the coordinate frame where k\\z and (61,82) — {eg,e^) tensor perturbations can be decomposed as ||l^ , [l^ , 

A^P{T,n,k) = [(1 - ^2)e2»^^i(fc) + (1 - ^,')e-'■^^e{k)] AP{T,f,,k) 
(A[f) +zA[J))(r,n,fc) = [(1 - fi)' e''^'^ e (k) + {1 + pi)^e-^'H\k)] A^p{T,^i,k) 
(Ag') - tA^p){T,n,k) = [{l + fi)'e'^'^e{k) + {1 - f^f e-^'H' {k)] A^^^(T,M,fc), 



(21) 



(22) 



where A^, and Ap are the variables introduced by Polnarev to describe the temperature and polarization pertur- 
bations generated by gravity waves. They satisfy the following Boltzmann equation @Jl^ 
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Aj, ' +ikfiA!^^ = -h - k[A^' > - ^ 



(T) 



(T) 



1 r 



(T) 



Ia 



(T) 



10 ™ ^ 7^T2 -r ^g-T4 5-PO -r ^-p2 



(23) 



Just like in the scalar case these equations can be integrated along the line of sight to give 

Af.^'Cro, n, k) = [(1 - ^i^)e^'^e{k) + (1 - ^?)e~^^H\k)\ T dre"'^4^^(fc, r) 

Jo 

(Ag^) + zA[J))(ro, n, fc) = [(1 - ^ife'^f'e{k) + (1 + A*)'e-2'*^2(^)] H dre^^f^ S^P {k, r) 

Jo 



(24) 



where 



s!^\k,T) = -he-'' 
S'-p{k,T) = -g^. 

Acting twice with the spin raising and lowering operators on the terms with gives 



(25) 



6\AP+^AP){ro,n,k)=e{k) 



'To 



dTsP\k,T) ( + 



1 - 



= £,\k)e^"^ / dTSP{k,T)[-£{x) - iB{x)] [(1 - ^?)e''=^'] 
Jo 

6 2(a(J) - iA^p){To,h, k) = e{k)e^'* £ drSPik, r) {^-df, - [(1 - /?)(! + y:fe'^^] 

= £}e^"^{k) / dTSP{k,T)[-£{x)+iB{x)][{\- ^l^)e''■'^'] 



(26) 

where we introduced operators £{x) = —12 + — 9^] — ^xdx and B{x) = 8a; + 2x'^dx- Expressions for the terms 
proportional to can be obtained analogously. 

('T'\ ('T'\ (T'^ 

For tensor modes all three quantities Aj. , A ~ and A ~ are non- vanishing and given by 



APiTo,n,k)= {I- fi')e'''^e{k) + il- fiyr'^feik) I drS).' > (t, k) e'' 
L J Jo 



To 



A^/^(ro,n,fc)= (l-/i2)g2»0^i(fc) + (i_^2)g-2.0^2(^) / 







Aj)(ro,n,fc)= [(l-/i^)e^'*fi(fc)-(l-M')e-^^^C'(fc)JS(x) dr^}.^ ^ (r, A;) e"^ 



To 



(27) 



From these expressions and equations (|^) , ( pl| ) one can explicitly show that B does not cross correlate with either T 
or E. 

The temperature power spectrum can be obtained easily in this formulation, 



^(T) _ 47r 



-J k^dkPh{k)Y^ J dnYi*^{n) " drSPik.T) {1 ~ ^i^)e^"^e'''^' 
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47r2lLl|ll / k^dkPhik) 



(1 + 2)1 



dTSP{k,T) / dfiP^i^i) (1 - fi^)e"'' 
J -I 



dTSP{k,T) I d^-^P,(^) (1-^2)2^^.;. 
-1 



477' 



= 47r 



{1 + 2)1 

2^-2)! 
I 



(An) 



(1 + 2) 

■ XI + 2)1 
{l^2)\ 



k^dkPhik) 
k^dkPhik) 
k^dkPh{k) 



dTSP{k,T) I d^ij-^Pi{^x){l + dlfe^^^ 
diiPii^i) {l + dl)\x^e^^^) 



dTS^p{k,T) 



/•To 

/ dTS!^f\k,T 

Jo 



(28) 



where we used Y,™ = [(2/ + \){l - m)\/{ATi){l + TO)!]i/2p™(^)e™'^ and PJ^i^i) = (-1)"(1 - fi^)""/^ £^Pi{n). Note 
that the calculation involved in the last step is the same as for the scalar polarization. The final expression agrees 
with the expression given in pl| ] , which was obtained using the radial decomposition of the tensor eigenfunctions [ p5| . 
Although the final result is not new, the simplicity of the derivation presented here demonstrates the utility of this 
approach and will in fact be used to derive tensor polarization power spectra. 

The expressions for the E and B power spectra are now easy to derive by noting that the angular dependence for 
A^-* and A^'' in (|2^) are equal to those for A^'' . The expressions only differ in the £ and B operators that can be 
applied after the angular integrals are done. This way we obtain using equation (|l 

2 



(Att)' 



k^dkPhik) 
k^dkPh{k) 
k^dkPhik) 







drSg\k^r)£{x)':^ 



drSg\Kr)\^m+j['ix)+^Jl^+'jM 



^ drSg\k^r)B{x)'M. 



= (4^)2 / ^dkPhik) ( / dTSg\k,T)\2j'i{x) + 



(29) 



For computational purposes it is convenient to further simplify these expressions by integrating by parts the deriva- 
tives j'i{x) and ji{x). This finally leads to 



A 



(T) 
E.Bl 



r drSg\k,r)^J^ 



d'rSgi{k,T)ji{x) 



sg\k,T)^g{^ 



kx 



2* 4^- 

A;^ kx 



2.972 



Sg\k,T)=, 



4^f 2* \ * 



(30) 



The power spectra are given by 



C 



(T) 
XI 



(47r)2 / k^dkPhik) A^,'(fc) 



eg) ^ (4vr)2 / edkP^{k)Ag\k)Ag\k), 



where X stands for T, E or B. Equations (|3C| ) and ( pi]) are the main results of this section 



(31) 



V. SMALL SCALE LIMIT 



In this section we derive the expressions for polarization in the small scale limit. The purpose of this section is to 
make a connection with previous work on this subject ^,^,^,0 and to provide an estimate on the validity of the 
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small scale approximation. In the small scale limit one considers only directions in the sky h which are close to 1, in 
which case instead of spherical decomposition one may use a plane wave expansion. For temperature anisotropics we 
replace 



so that 



To expand s — ±2 weighted functions we use 

■{1-2)1 ' 



iie 



(32) 



(33) 



-2Ylr 



{1 + 2)1 
(1-2)1 



{I 



{3'n 



(2^)-2l^"2gJ.0^ 



(34) 



which leads to the following expression 

(Q + iU){n) = -(2^)2 J dH [E{1) + iB{l)]^d ^e'^ ^ 
{Q - iU){n) = -(2^)2 J dH [E{1) ~ iB{l)]^6 ^e'^'^. 



(35) 



From equation (A2) wc obtain in the small scale limit 



-2i{,p-4>i) d-O 



(36) 



where {l^ + Hy) = ^e*"^'. 

The above expression was derived in the spherical basis where ei = eg and 62 = e^, but in the small scale limit 
one can define a fixed basis in the sky perpendicular to z, e'l = e.^ and 63 = e^. The Stokes parameters in the two 
coordinate systems are related by 



{Q + iUy = e-^'^{Q + iU) 
{Q-iUy = e^'^{Q-iU). 



(37) 



Combining equations 



57|) we find 
Q'{e) = (27r)-2 



d^Z [E{1) cos{2(f>i) - B{1) iiin{2cj)i)]e 



le 



U' {9) = {2n)-'^ J d^l [E{l)sm{2(t>i) + B{l)cos{2<j)i)]e'''-^. 



(38) 



These relations agree with those given in |14|, which were derived in the small scale approximation. As already shown 
there, power spectra and correlation functions for Q and U used in previous work on this subject can be 

simply derived from these expressions. Of course, for scalar modes B'-^^I) — 0, while for the tensor modes both 
E^'^'>{1) and B^^\l) combinations contribute. 

The expressions for Q and U (equation 38 ) are easier to compute in the small scale limit than the general expressions 
presented in this paper (equation because Fourier analysis allows one to use Fast Fourier Transform techniques. In 
addition, the characteristic signature of scalar polarization is simple to understand in this limit and can in principle 
be directly observed with the interferometer measurements p4| . On the other hand, the exact power spectra derived 
in this paper (equations |l^ and |3^, ^) are as simple or even simpler to compute with the integral approach than 
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their small scale analogs. Note that this need not be the case if one uses the standard approach where Boltzmann 
equation is first expanded in a hierarchical system of coupled differential equations [Q . In Fig. ^ we compare the exact 
power spectrum (solid lines) with the one derived in the small scale approximation (dashed lines), both for scalar E 
(a) and tensor E (b) and B (c) combinations. The two models are standard CDM with and without reionization. The 
latter boosts the amplitude of polarization on large scales. The integral solution for scalar polarization in the small 
scale approximation was given in and is actually more complicated that the exact expression presented in this 
paper. In the reionized case the small scale approximation agrees well with exact calculation even at very large scales, 
while in the standard recombination scenario there are significant differences for I < 30. Even though the relative 
error is large in this case, the overall amplitude on these scales is probably too small to be observed. 

For tensors the small scale approximation results in equation ( ^0| ) without the terms that contain or x~^. 
Because ji{x) ~ for a; < Z these terms are suppressed by l^^ and l^^, respectively, and are neghgible compared 
to other terms for large I. The small scale approximation agrees well with the exact calculation for B combination 
(Fig. |l|c), specially for the no-reionization model. For the E combination the agreement is worse and there are 
notable discrepancies between the two even at I ^ 100. We conclude that although the small scale expressions for 
the power spectrum can provide a good approximation in certain models, there is no reason to use these instead 
of the exact expressions. The exact integral solution for the power spectrum requires no additional computational 
expense compared to the small scale approximation and it should be used whenever accurate theoretical predictions 
are required. 



VI. ANALYSIS OF ALL-SKY MAPS 



In this section we discuss issues related to simulating and analyzing all-sky polarization and temperature maps. 
This should be specially useful for future satellite missions which will measure temperature anisotropics and 

polarization over the whole sky with a high angular resolution. Such an all-sky analysis will be of particular importance 
if reionization and tensor fluctuations are important, in which case polarization will give useful information on large 
angular scales, where Fourier analysis (i.e. division of the sky into locally flat patches) is not possible. In addition, it 
is important to know how to simulate an all-sky map which preserves proper correlations between neighboring patches 
of the sky and with which small scale analysis can be tested for possible biases. 

To make an all-sky map we need to generate the multipole moments or^/m, aE^im and aB i„^. This can be done by a 
generalization of the method given in jl^. For each I one diagonalizes the correlation matrix Mn = Cti, M22 — Csh 
M12 = M21 = Cci and generates from a normalized gaussian distribution two pairs of random numbers (for real and 
imaginary components of ai±rn.)- Each pair is multiplied with the square root of eigenvalues of M and rotated back 
to the original frame. This gives a realization of aT.i±m and aE,i±m with correct cross-correlation properties. For 
o-B,i±m the procedure is simpler, because it does not cross-correlate with either T or i5, so a pair of gaussian random 

1/2 

variables is multiplied with C^j to make a realization of aB,i±m- Of course, for scalars aB.im = 0. 

Once aE,im and aB,im are generated we can form their linear combinations a2,im and a_2,;m, which are equal in 
the scalar case. Finally, to make a map of Q{fi) and U{fi) in th e sky we perform the sum in equation (|^), using 



the explicit form of spin- weighted harmonics sYim{'f>') (equation A6). To reconstruct the polarization power spectrum 
from a map of Q{n) and C/(n) one first combines them i\iQ + iU and Q — ill to obtain spin ±2 quantities. Performing 
the integral over ±2Yim (equation ||) projects out ±20/™, from which aE,im and aB,im can be obtained. 

Once we have the multipole moments we can construct various power spectrum estimators and analyze their 



variances. In the case of full sky coverage one may generalize the approach in |26| to estimate the variance in the 
power spectrum estimator in the presence of noise. We will assume that we are given a map of temperature and 
polarization with Npi^ pixels and that the noise is uncorrelated from pixel to pixel and also between T, Q and U. 
The rms noise in the temperature is ctt and that in Q and U is ap. If temperature and polarization are obtained 
from the same experiment by adding and subtracting the intensities between two orthogonal polarizations then the 
rms noise in temperature and polarization are related by (t|, — crp/2 [Q. 

Under these conditions and using the orthogonality of the sYim we obtain the statistical property of noise, 

//'^noisoN*^noisc \ _ 477(7 j 

//'_noise\*„noisc \ OTTCTp 

\(«2,ImJ "2, Cm'/ - -ITj OU'<^mm' 

pix 

Q 2 

//^noisc \* noise \ " P x X 

(l«-2,;mJ a-2,;'m'/ " ^T^ ^U'f^vara' 
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//■^noise \* noise 



0, 



(39) 



where by assumption there arc no correlations between the noise in temperature and polarization. With these and 
equations (|6|j|) we find 





= {Crie- 


2^2 




{'^*E,lm'^E,l'm') 


= [Cme- 






{(^*B,lmO'B,l'm') 


= {Csie- 




' Wp^)Sll'Sjnrn 






''Oil' 




{0'*B,l'm'0'E,lm) 




aT,ln 


>=0. 



(40) 

For simplicity we characterized the beam smearing by e' '^''/^ where at is the gaussian size of the beam and we defined 
The estimator for the temperature power spectrum is pq], 



G 



Tl 



E 



21 + 1 



(41) 



Similarly for polarization and cross correlation the optimal estimators are given by [p^ 



C 



El 



c. 



Bl 



O-BJ-n 



Cci — 



rn 

E 

m 

E 



Wr 



21 + 1 

2 



Wr 



21 + 1 

{a*E,hn^T.lm. + aEJmO^, ;„J 



2{2l + l) 



(42) 



The covariance matrix between the different estimators, Cov{XX') = {{X — {X)){X' — {X'))) is easily calculated 
using equation (|o|). The diagonal terms are given by 

[CTl+W^^'^^'^f 

{CEi+Wp^e'"''"f 
{CBi+Wp^e''<f 



Cov((7|,) = 
Cov((7|;) = 
Cov((7^;) = 

The non-zero off diagonal terms are 



21 + 


1 


2 




21 + 


1 


2 




21 + 


1 


1 




21 + 


1 



Cli + {Cti + u;^^e''"')(C'B/ + u-p^e''"') 



(43) 



Co^{CtiCei) 
CoNiCriCci) 



2Z + 1 
2 



^^^-^Cci(CTi+w-T,'e' 
Cov{CeiCci) = ^j^'^ciiCEi + wp^e''^'). 



(44) 



These expressions agree in the small scale limit with those given in |T^]. Note that the theoretical analysis is more 
complicated if all four power spectrum estimators are used to deduce the underlying cosmological model. For example, 
to test the sensitivity of the spectrum on the underlying parameter one uses the Fisher information matrix approach 
. If only temperature information is given then for each I a derivative of the temperature spectrum with respect to 
the parameter under investigation is computed and this information is then summed over all I weig htedbyCov-^Cf,,). 
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In the more general case discussed here instead of a single derivative we have a vector of four derivatives and the 
weighting is given by the inverse of the covariance matrix, 



"^.-EE^Cov-(^-^-)^, (45) 

I X,Y 

where is the Fisher information or curvature matrix, Cov~^ is the inverse of the covariance matrix, Si are the 
cosmological parameters one would like to estimate and X, Y stands for T, i?, B, C. For each I one has to invert the 
covariance matrix and sum over X and Y , which makes the numerical evaluation of this expression somewhat more 
involved. 



VII. CONCLUSIONS 



In this paper we developed the formalism for an all-sky analysis of polarization using the theory of spin-weighted 
functions. We show that one can define rotationally invariant electric and magnetic- type parity fields E and B 
from the usual Q and U Stokes parameters. A complete statistical characterization of CMB anisotropics requires 
four correlation functions, the auto-correlations of T, E and B and the cross-correlation between E and T. The 
pseudo-scalar nature of B makes its cross-correlation with T and E vanish. For scalar modes B field vanishes. 

Intuitive understanding of these results can be obtained by considering polarization created by each plane wave 
given by direction k. Photon propagation can be described by scattering through a plane-parallel medium. The 
cross-section only depends on the angle between photon direction n and k, so for a local coordinate system oriented 
in this direction only Q Stokes parameter will be generated, while U will vanish by symmetry arguments | p2[ |. In the 
real universe one has to consider a superposition of plane waves so this property does not hold in real space. However, 
by performing the analog of a plane wave expansion on the sphere this property becomes valid again and leads to the 
vanishing of B in the scalar case. For tensor perturbations this is not true even in this fc dependent frame, because 
each plane wave consists of two different independent "polarization" states, which depend not only on the direction 
of plane wave, but also on the azimuthal angle perpendicular to fc. The symmetry above is thus explicitly broken. 
Both Q and U are generated in this frame and, equivalently, both E and B are generated in general. 

Combining the formalism of spin-weighted functions and the line of sight solution of the Boltzmann equation we 
obtained the exact expressions for the power spectra both for scalar and tensor modes. We present their numerical 
evaluations for a representative set of models. A numerical implementation of the solution is publicly available and 
can be obtained from the authors |2^]. We also compared the exact solutions to their analogs in the small scale 
approximation obtained previously. While the latter are accurate for all but the largest angular scales, the simple 
form of the exact solution suggests that the small scale approximation should be replaced with the exact solution for 
all calculations. If both scalars and tensors are contributing to a particular combination then the power spectrum for 
that combination is obtained by adding the individual contributions. Cross-correlation terms between different types 
of perturbations vanish after the integration over azimuthal angle (j) both for the temperature and for the E and B 
polarization, as can be seen from equations (15) and (p7|). This result holds even for the defect models, where the 
same source generates scalar, vector and tensor perturbations. 

In summary, future CMB satellite missions will produce all-sky maps of polarization and these maps will have to be 
analyzed using techniques similar to the one presented in this paper. Polarization measurements have the sensitivity 
to certain cosmological parameters which is not achievable from the temperature measurements alone. This sensitivity 
is particularly important on large angular scales, where previously used approximations break down and have to be 
replaced with the exact expressions for the polarization power spectra presented in this paper. 
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APPENDIX A: SPIN- WEIGHTED FUNCTIONS 



In this Appendix we review the theory of spin- weighted functions and their expansion in spin-s spherical harmonics. 
This was used in the main text to make an all-sky expansion of Stokes Q and U Stokes parameters. The main 
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Our 

We 



application of these functions in the past was in the theory of gravitational wave radiation (see e.g 
discussion follows closely that of Goldberg et al. jl^, which is based on the work by Newman and Penrose 
refer to these references for a more detailed discussion. 

For any given direction on the sphere specified by the angles {0,(f>), one can define three orthogonal vectors, one 
radial and two tangential to the sphere. Let us denote the radial direction vector with n and the tangential with ei, 
62. The latter two are only defined up to a rotation around n. 

A function sf{d^(j)) defined on the sphere is said to have spin-s if under a right-handed rotation of (61,62) by 
an angle ■0 it transforms as s/'(0, </>) = e~'^^^ sf{d,(f>)- For example, given an arbitrary vector a on the sphere the 
quantities a • ei + ia • 62, n • a and a • ei — ia • 62 have spin 1,0 and -1 respectively. Note that we use a different 
convention for rotation than Goldberg et al. ||l9[| to agree with the previous literature on polarization. 

A scalar field on the sphere can be expanded in spherical harmonics, Yi^id, (p), which form a complete and orthonor- 
mal basis. These functions are not appropriate to expand spin weighted functions with s ^ 0. There exist analog sets 
of functions that can be used to expand spin-s functions, the so called spin-s spherical harmonics sYimid, </>). These 
sets of functions (one set for each particular spin) satisfy the same completness and orthogonality relations. 



sY;^{d,cj)) sYiUO' ,cl)') ^ 6{cl> - cl>')6{cose - cose'). 



(Al) 



Im 



An important property of spin-s functions is that there exists a spin raising (lowering) operator d {d ) with the 
property of raising (lowering) the spin-weight of a function, {d sfY = e^*(*+i)'/'(5 j,/, ^f)' — Their 
explicit expression is given by 



9 rn^9 



sin-^(0) ,/(0,0) 



(A2) 



In this paper we are interested in polarization, which is a quantity of spin ±2. The d and d operators acting twice 
on a function ±2f{fJ-, 0) that satisfies sf = ifn sf can be expressed as 



2/(^,0) 



l-/i^ 



[(1-^*2)2/(^,0)] 

[(l-A*2)_2/(^,0)] 



(A3) 



where /i = cos{6). With the aid of these operators one can express sYim in terms of the spin zero spherical harmonics 
Yira, which are the usual spherical harmonics. 



il + s)\ 

(i + sy. 



&'Yirn ,(0<S<0 

{^iyd-'Yi„, ,{-i<s<o). 



(A4) 



.(/-.)!_ 

The following properties of spin-weighted harmonics are also useful 

s^im — {—ly -sYl-,n 
d sYlra = [{1 - -syi + S + 1)]^ ,+ iYim 
d sYlr,^ =-[(? + m - S + 1)]* 

3 d sYi„, ^-{l-s){l + s + l) sYi„^ (A5) 
Finally, to construct a map of polarization one needs an explicit expression for the spin weighted functions. 
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{I + m)\[l - m)\ 21 + I 



E 



r + s 



{l + s)\{l~s)\ 47r 



cot 



1/2 

sin2'(e'/2) 



(A6) 
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FIG. 1. Comparison between exact calculation (solid lines) and small scale approximation (dashed lines) for standard CDM 
model with and without reionization. In the latter case we use optical depth of 0.2. The reionized models are the upper curves 
on large scales. The comparison is for scalar E (a) and tensor E (b) and B (c) polarization power spectra. The spectra are in 
units of Tq — {2.729K)^ and are normalized to COBE. While the predictions agree for large I there are significant discrepancies 
in certain models for small I. 
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